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In this study, we present the effects of functionahzation of single carbon nanotubes (CNT) with 
simple organic molecules -NH, -NH2, -CH2, -CH3, -OH, -COOH and both -COOH, -OH, 
analyzing elastic properties of the structures. We have performed DFT calculations for CNT with 
various densities of the attached molecules. We have determined the changes in the geometry, the 
Young's modulus. Shear modulus and Bulk modulus and Poisson's ratio as a function of the density 
of the adsorbed molecules. For better comparing with some experimental works we have calculated 
percentage change of those quantities due to pristine CNT. Increasing number of covalently bound 
groups to the lateral surface of CNT decrease Young's modulus. Shear modulus and Bulk modulus. 
The biggest relative diffrences for each of those quantities, for 12,5% concentration, are ca. 30% for 
structures functionalized with —CH2. For —NH it is similar, in contrast to Shear modulus where 
reduction is smaller and around 15 %. For remaining structures, we have noted 20 % reduction in 
Young's and Bulk modulus, and 10 % in Shear modulus. Whereas Poisson's ratio is little changed 
in order to pristine CNTs. 

PACS numbers: 31.15.A-, 31.15.ae, 31.15.E-, 31.15.ec, 61.46.-w, 61.46.Fg, 61.48.De, 81.05.U-, 81.07.De 



I. INTRODUCTION 

Remarkable electronic, mechanical and thermal prop- 
erties of CNTs have been developed rapidly since CNT 
was discovered in 1991 by lijima^. New composite mate- 
rials made by adding CNT to various materials such as al- 
loys, polymers, and metals are expected to enhance their 
properties, such as mechanical strength, electrical and 
thermal conductivity, and chemical stability^^^. How- 
ever, pristine CNTs are not soluble in water or in organic 
solvents and have tendencies to create bundles. This 
drastically limits its use in industrial applications like 
nano-composites. Only simple organic molecules (such 
as CHn, NHn fragments, OH, COOH groups) adsorbed 
at the surface of CNTs allow for strong binding of the 
functionalized in such a way CNTs with matrix material, 
typically a polymer or a metal^^. This is the possible 
answer for problems with dispersion of CNTs and their 
interaction with matrix. 

From the other point of view, functionahzation to the 
side wall of CNTs introduce changes in structure, and 
defects^^^, which could decrease the strenght of the 
structure due to pristine CNT. Having into considera- 
tion other possibly applications for functionalized CNTs, 
for instance in microeletronics or medicine, it is crucial to 
anlyzed elastic properties of those structures. Compar- 
ing the extensive amount of works of elastic properties of 
pristine CNTs with few papers concerning functionalized 
CNTs there is a need to study and predict the influence 
of functionahzation on CNTs. 

In this study, we have analized the effects of func- 
tionahzation of single wall (9,0), (10,0) and (11,0) CNTs 
(SWNTs) with simple organic fragments CHn, NHn, OH, 
and COOH, focusing on their Young's, Shear and Bulk 
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FIG. 1. Example structures of funtionalized CNT with simple 
(-NH, -NH2 , -CH2, -CH3, -COOH, -OH, -OH and -COOH) 
organic fragments after fully optimization. 



modulus. 



II. CALCULATION DETAILS 



We have analyzed elastic properties of three types of 
CNT: (9,0), (10,0) and (11,0). All of CNT have been 
covalently functionalized, up to 12.5 % concentration, 
by attaching to their lateral surface simple organic frag- 
ments such as -NH, -NH2, -CH2, -CHs, -OH, 
-COOH and both -COOH, -OH (Fig.H]). 
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The total energies and components of stress tensor are 
obtained from the ab initio calculations in the framework 
of the density functional theory^^^^^. Stability studies in- 
cluding analysis of binding and adsorption energies due 
to number of attachments and radius of CNTs were pre- 
sented in our previous published works^"— . Therefore 
here we have decided to show only adsorption energy per 
molecule, called also packing energy by Cotc^^. 
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We use the generalized gradient approximation (GGA) 
for the exchange-correlation density functional^^ and 
supercell geometry within the numerical package 
SIESTAi^ii^. A kinetic energy cut-off (parameter Mesh- 
Cutoff in the SIESTA code) of 300 Ry and split double 
zeta basis set with spin polarization have been used in 
all calculations. Each supercell contains two primitive 
unit cells. The lateral separation (i.e., lateral lattice con- 
stants in the direction perpendicular to the symmetry 
axis) has been set to 30 A, just to eliminate completely 
the spurious interaction between neighboring cells. 

The lattice constant along the symmetry axis of the 
tube has been also optimized. The positions of all atoms 
in the supercell have been optimized to get vanishing 
forces on atoms. This provides us also with the Poisson 
ratio values, which has been calculated using following 
equation: 



Ar / 



(2) 



where Ar is change of radius in strained nanotube. 

The most interesting quantity. Young's modulus, have 
been determined in two ways: (i) - by comparing the total 
energy of unstressed {Ei) and stressed (Ei-^^i) systems, 
applying the tensile stress along the symmetry axis of 
CNTs: 
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where / is a lattice constant along the axis of the tube, 
Al is elongation in the chosen direction, and Vq is volume 
of unstressed system; and (ii) - from components {(Ju) of 
the stress tensor: 



Y = ^ 



(4) 



Volume of the pure CNT has been calculated using fol- 
lowing equation: 



Vn = 2-7r -r -I't, 



(5) 



where thickness t has been chosen as double Van der 
Waals radius of C atom (0.34 nm)^^^^ and r being nan- 



otube's radius calculated as geometrical average of car- 
bon atom positions on the CNT surface. In case of func- 
tionalized CNT we have neglected volume of attached 
molecules. 

Bulk modulus (Eq. [6]) and Shear modulus (Eq. [7j) have 
been also calculated according to the formulas: 
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III. RESULTS AND DISCUSSION 
A. Influence of functionalization on the structure 

In previous works^^^ we have shown how functional- 
ization with simple molecules influences structural and 
electronic properties of CNT. Therefore, in this paper, 
we present only the strenght of covalent binding beetwen 
CNT and its attachements, which is described in term 
of adsorption energy per molecule. This is the amount 
of energy which is necessery for attachment of each re- 
garded fragment to lateral surface of CNT. Adsorption 
eneregy per molecule, which is shown in Fig. [21 allows to 
distinguish between radicals (-NH and -CH2) and non- 
radicals (-NH2, -CH3, -COOH and -OH,). Ah of the 
fragments induce rehybridization from sp^ to sp^ of C-C 
bonds in neighbourhood of attachemnt, but some of them 
also cause other deformation of structure. Stronger co- 
valent binding which means lower adsorption energy per 
molecule leads to more changes in structure. 

We have compared in Fig. [3] two types of frag- 
ments: -OH groups nad -CH2 radicals. On the contrary 
to hydroxyl group, which leave CNT structure almost 
uchanged, -CH2 for some concentration could destroy 
whole CNT. Adsorption energy per molecule for radicals, 
like -CH2, which are stronger bound to the lateral surface 
of CNT then -OH groups are almost indepedent of diam- 
eter of the tube. This quantity, for radicals, is slightly 
increasing with increasing density of attached molecules. 
For other non-radical attachments adsorption energy per 
molecule is maintained on the same level for all of re- 
garded concentrations. 

Functionalization has not only local influence on struc- 
ture of CNT, but is also global. Radius of CNT rapidly 
changes even for small concentration of radicals attached 
to the lateral surface (Fig. |4]). For example, radius of 
pure (9,0) CNT equals 3.592 A and for CNT with 9 -CH2 
- 3.715 A, which means 3.31% change. For other groups 
we have also observed an increasing trend in comparision 
to the pure case. Percentage change for -CH3 is 0.92 %, 
approximately 3.6 times smaller than for -CH2. Radius 
of CNT, likewise adsorption energy per molecule, varies 
with the strenght of binding. 
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FIG. 2. (color online) Adsorption energy per molecule of (9,0) 
CNT functionalized with -NH, -NHs , -CH2, -CH3, -COOH, 
-OH, -OH and -COOH groups as a function of number of 
covalently bound fragments to sidewall of tube. Number of 
groups (bottom axis) have been recalculated to percentage 
concentration (top axis). 
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FIG. 3. (color online) Adsorption energy per molecule of 
(9,0), (10,0), (11,0) CNT functionalized with -CH2 and -OH 
fragments as a function of density of attached molecules (a) 
and percentage concentration (b) . Number of groups (bottom 
axis in (b)) have been recalculated to percentage concentra- 
tion (top axis in (a)). 



FIG. 4. (color online) Radius of (9,0) CNT functionalized 
with -NH, -NH2 , -CH2, -CH3, -COOH, -OH, -OH and - 
COOH groups as a function of number of covalently bound 
fragments to sidewall of tube. Number of groups (bottom 
axis) have been recalculated to percentage concentration (top 
axis). Radius of functionalized CNT is increasing with in- 
creasing number of attached fragments. 
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FIG. 5. (color online) Radius of (9,0), (10,0) and (11,0) CNT 
functionalized with -CH2 and -OH as a function of number of 
attached fragments. Top axis with percentage concentration 
refers to (11,0) CNT, middle - (10,0) CNT, and bottom - (9,0) 
CNT. 



As mentioned before, attaching -CH2 radicals to the 
lateral surface of CNT causes strong deformations like 
5-7 defects. Functionalization with -OH groups slightly 
changes the cross-section of CNT - from circle to ellipse. 
Therefore we have decided to compare (Fig. [5]) both types 
of attachments for all of the CNT - (9,0), (10,0) and 
(11,0). We have noticed, for biggest considered concen- 
tration - 12.5%, that (9,0) CNT functionalized with -CH2 
shows the biggest percentage change ( 3.31%) of radius 
with comparision to (10,0) (2.83%) and (11,0) (1.81%) 



CNT. For -OH groups the same trend, but weaker per- 
centage decrease is observed. The percentage changes 
are 0.64%, 0.43% and 0.39%, for (9,0), (10,0) and (11,0) 
CNT, respectively. Therefore, one can deduce that func- 
tionalization of bigger nanotubes has less influence on its 
structure than functionalization CNT with smaller diam- 
eter. 
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B. Elastic properties of pure CNT 

For the purpose of comparison our calculations to pre- 
vious works and to get the point of reference to function- 
alized case we have analyzed elastic properties of pure 
(9,0), (10,0) and (11,0) CNTs. 

TABLE L Elastic properties of (9,0), (10,0) and (11,0) CNTs. 



Property 


(9,0) 


(10,0) 


(11,0) 


Y (TPa) 


1.02 


1.03 


1.02 


K (TPa) 


0.61 


0.57 


0.54 


G (TPa) 


0.41 


0.43 


0.43 




0.22 


0.20 


0.18 



Young's moduli. Shear moduli. Bulk moduli, and Pois- 
son's ratios for (9,0), (10,0) and (11,0) CNTs are gathered 
in TablH The calculated Young's moduli of the pure CNT 
compare excellently to experimental findings (0.32 TPa - 
1.80 TPa)^i^i^ and previous theoretical works (0.8 - 1.5 
TPa)^ii^i^i^i^"— . Calculated Poisson ratios are iden- 
tical to the experimental ones and very close to theo- 
retical predictions (0.19-0.34)^i^i^i^, and theoretical 
and experimental values of the Shear modulus and Bulk 
modulus are 0.45 TPa-0.58 TPa^^^^^^i 23,30 q 5Q r^p^^ 

-0.78 TPaP^iii, respectively. 

We have also checked dependence of all the elastic 
properties for wider range of zigzag CNTs. Only for 
small CNT, like (4,0) and (5,0) all the values are smaller. 
For bigger CNT, up to (20,0), the values are very sim- 
ilar tthose shown in TabUl From (6,0) CNT ah of the 
elastic moduli become insensitive to diameter of CNT. 
This behaviour has been remarked for Young's modulus 
b3^^^^-3i fo^ ^Yie^^ modulus byl2^^. 



C. Elastic properties of functionalized CNT 

In this part we have performed study of elestic prop- 
erties of functionalized CNT. 

First of all we have analyzed Young's modulus of (9,0) 
CNT functionalized with -NH, -NH2 , -CH2, -CH3, - 
COOH, -OH, -OH and -COOH groups (Fig. For ah 
of groups this quantity is decreasing with increasing den- 
sity of attachments. However, for radicals, -CH2 and 
-NH, trend is much more stronger then for other groups. 
CNT functionalized with 9 CH2 (12.5 %) has reduced 
its structure by 28.41 %. Whereas CNT with 9 -CH3 
groups reduction in Young's modulus is equal 13.52 %. 
Changes in structure cause by functionalization, which 
have been just described in terms of CNT radius (Fig.|4]), 
are also very well noticable in Young's modulus. Consid- 
ering that similar dependence for armchair CNTs func- 
tionalized with -COOH was reported by Cotc^^ this trend 
can be generalized for all of SWCNT. 

For the purpose of comparing how Young's modulus 
depends of diameter of the tubes we have chosen -OH 
groups and CH2 fragments. In Fig. [3 we have plotted the 
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FIG. 6. (color online) The Young's modulus of (9,0) CNT 
functionalized with -NH, -NH2 , -CH2, -CH3, -COOH, -OH, 
-OH and -COOH groups as a function of density of attached 
molecules. Number of fragments bound to CNT have been 
recalculated to percentage concentration (top axis). On the 
right axis we have depicted percentage change of Young's 
modulus due to increasing functionalization. Young's mod- 
ulus is decreasing with increasing number of attached groups 
(concentration). 
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FIG. 7. (color online) The Young's modulus of (9,0), (10,0) 
and (11,0) CNT functionalized with -OH , and -CH2 frag- 
ments as a function of density of attached molecules (a) and 
percentage concentration (b). Young's modulus is almost in- 
dependent of diameter of the tube. 



dependence Young's modulus of (9,0), (10,0) and (11,0) 
CNT functionalized with -OH and -CH2 on density . - 
OH groups represents behaviour of non-radical attache- 
ments which cause small deformation of CNT, and -CH2 
- radicals wich induce 5-7 defects and for bigger concen- 
trations are able to destroy whole structure. It is clearly 
seen (Fig. [7^) that for -CH2 Young's modulus is slightly 
dependent and for -OH we can not distinguish between 
tubes. 
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FIG. 8. (color online) The Shear modulus of (9,0) CNT func- 
tionalized with -NH, -NH2 , -CH2, -CH3, -COOH, -OH, - 
OH and -COOH groups as a function of density of attached 
molecules. Number of fragments bound to CNT have been 
recalculated to percentage concentration (top axis). On the 
right axis we have depicted percentage change of Shear mod- 
ulus due to increasing functionalization. Shear modulus, like 
Young's modulus is decreasing with increasing number of at- 
tached groups (concentration). 



FIG. 9. (color online) The Bulk modulus of (9,0) CNT func- 
tionalized with -NH, -NH2 , -CH2, -CH3, -COOH, -OH, - 
OH and -COOH groups as a function of density of attached 
molecules. Number of fragments bound to CNT have been 
recalculated to percentage concentration (top axis). On the 
right axis we have depicted percentage change of Young's 
modulus due to increasing functionalization. Bulk modulus, 
like Bulk modulus and Shear modulus is decreasing with in- 
creasing number of attached groups (concentration). 



Poisson ratio for structures functionalized by all of the 
fragments always oscilate between 0.17 and 0.24. This 
quantity, for this range of concentration, do not allow 
to distinguish among different types of attachments and 
CNTs. Hovewer Poisson ratio is necessity for calculating 
other elestic moduli. 

Shear modulus has been plotted as a function of num- 
ber and concentration of attached molecules (Fig. [8]). 
Shear modulus drops slightly further, if one increases 
the density of the attached molecules, especially non- 
radicals, to the CNTs. For the -COOH groups in con- 
centration less then 5.56 %(4 groups) we have noted ex- 
cellent agreement with Zheng^^, shear modulus is very 
little change (less then 4.63 %). 

Only for non-radical attachments in less then 10 % 
concentration, we have noted reduction in Shear modulus 
little higher then 10%. We have proved that Franklad's^^ 
prediction of negligible influence functionalization on 
shear modulus can not be generalized to all of simple 
attachemnets. 

The last elastic moduli, which we have taken into con- 
sideration, is Bulk modulus. This modulus as function 
of number of attached molecules and pecentage concen- 
tration, likewise Young's and Shear moduli, is shown in 
Fig. m Functionalize CNTs repsonse to uniform pres- 
sure likewise to linear strain and to shearing strains. The 
percentage reduction in Bulk modulus for non-radical at- 
tachemnets is greater then in Shear and Young's modu- 
lus. 



IV. CONCLUSIONS 

We have performed extensive studies of the elastic 
properties of the various types of CNTs functionalized, 
up to 12.5 % concentration, with many molecular groups. 
We have analized the crucial elastic properties such as 
Young's modulus, Poisson ratio, Shear modulus and Bulk 
modulus. 

We have shown by analyzing adsorption energy per 
molecule and radius of nanotube that functionalization 
to the lateral surface change the whole structure. Re hy- 
bridization of C-C bonds, deformation and defects have 
significant influence on elastic properties. Stronger co- 
valent bond leads to larger CNT deformation, resulting 
in decreasing of Young's, Shear and Bulk modulus. The 
Bulk modulus of functionalized CNTs, in contrast to - 
CH2 case, reveal the greatest percentage changes in com- 
parision to pure CNT. 

As far as Young's, Shear and Bulk modulus reflect 
changes in structure, the Poissons ratio remains almost 
unchanged. Those slight changes have no obvious ten- 
dency. 

CNT covalently functionalized with small concentra- 
tion non-radical groups could be reinforcement for com- 
posite materials. Functionalization with -COOH, -OH, 
-NH2, -CH3 do reduce all of the elastic moduli, but not 
to the extent that could disturb strong enforcement of 
the composite materials. The reduction is not too high - 
approximately 15 % in Young's modulus, 10 % in Shear 
modulus and 20% in Bulk modulus in comparision to 
pure CNTs. From the other hand they are indispensable 
to binding CNT to polymer matrix and significantly im- 



6 



prove homogenous dispersion and integration of CNTs in 
polymers, reduce its tendency to reagglomeration. This 
feature substantialy enhence elastic properties of poly- 
mer matrices with CNTs, even more then functionaliza- 
tion spoil because of inducing changes in tube's struc- 
ture resulting reduction in elastic moduli. This effect 
regarding to Young's modulus has been already exper- 
imentaly observed^^^^^^ for amines,^ amides, and car- 
boxyhc groups^i^. But for CNT with hydroxyl groups 
dispersed into polymer matrix Wang^ - reported reduc- 
tion in Young's modulus in comparition to pure matrix. 
Taking into account all of the possible applications, not 



only nano-composites, we can conclude that influence of 
functionalization on elastic properties of those structure 
can not be neglected. 
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